EG21007 - INTRODUCTION TO PROGRAMMING 
Assignment 4: Hanging elastic chain 
Due by: Wednesday, 19 th Nov 2014, 23:59 



A horizontally placed chain consists of N point particles of mass m, connected to each other 
by massless harmonic springs with the spring constant k. Let us assume that both ends of 
the chain are fixed and cannot move. If the chain is placed in the field of gravity it will deform 
under its own weight. The shape of the deformed chain will minimize the total energy (sum of 
elastic and gravitational energies). This can serve as a very crude model of a hanging 
bridge. 

In Figure 1 (upper panel), we show the shape of the undeformed chain (before gravity was 
switched on). In the lower panel we show shape of the deformed chain. Please note that 
first and the last bead are fixed and cannot move. 

If we assume that particles can only move in the vertical (y) direction and that the relative 
displacement between neighbouring beads is small (large spring constant), total energy of 
the system can be approximated as: 
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where g = 9.81 m/s 2 is the acceleration of gravity and we have omitted unimportant constant 
terms. 
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Figure 1: Chain of point-like particles of masses m connected with massless harmonic springs with spring 
constant k. Top: Undeformed chain. Bottom: Shape of the chain under its own weight (i.e., after the gravity was 

switched on). 

Task: 

Write a program that: 

a) finds optimal shape of a hanging chain by minimizing the energy function defined 
above; 

b) plots the optimal shape; 



i-i, 



c) calculates and prints on screen relative change in length, defined as S = 

/in the length of the deformed chain and b is the length of the initial (straight) 
configuration; 



where 



Assume that number of particles is N=1 01 , m,=0.5kg (500g) and /c=10 4 N/m 2 . Also assume 
that in the initial configuration particles are placed in a straight line with the 0.1m (10cm) 
spacing, i.e., Xi=0.1i (i=1,...,N). Label both axes appropriately. 

Remember, first and the last particles cannot move! Python indices start from 0, so be 
careful when implementing energy function into your code! 

Hint: You can use function fminbfgs in the scipy. optimize module to find the minimum 
energy configuration. Assume that the initial configuration (xO in the fmin bfgs notation) is a 
NumPy array on length N filled with zeros. 

NOTE: 

Half the available marks for this assignment will be awarded for originality, creativity, style, 
readability, presentation, quality of documentation, and design features incorporated into the 
code. Thus in general it will be insufficient to simply satisfy the requirements as outlined for 
this assignment with the minimum of thought and effort. Excellent code, however, will be 
awarded an excellent mark. Source code in Python 3 that satisfies the requirements of the 
assignment is to be submitted online by way of the my. Dundee web portal (within 
EG21007^ [Assignments]) by 11:59pm Wednesday 19 th November 2014. 

You should submit two files. One is the Python code and the other is a brief PDF document 
outlining what the program does, what parameters were used to test it and any potential 
problems you encountered. This should not be a full report, just a few sentences. However, it 
should be coherent, properly formatted and self-contained. 

Your Python solution to this assignment should be submitted as a file identified by your 
name: 

assignment_4_code_J_Smith.py 

where "J" would be replaced by your first initial (say, "E" for Emily) and "Smith" should be 
replaced by your last name. 

The PDF document for this assignment should be identified by your name: 

assignment_4_writeup_J_Smith.pdf 
with the same naming rules as above. 

A submission that does not comply with the prescribed naming convention will carry a 5% 
penalty. 

In accordance with the University policy on late submissions, assignments submitted after 
1 1 :59pm Wednesday 19 th November 2014 will attract an immediate 20% penalty. 
Assignments more than one week late will not be considered without a medical certificate 
that provides acceptable reason for the delay. 

Collaborative work: Note that collaborative work is not allowed. While you are welcome to 
discuss this assignment with each other, each students is required to write the program by 
him/herself. If I encounter programs that are clearly similar to each other (cosmetic changes 
such as renaming variables and/or modifying comments do not count as different programs!) 
all involved students will have their grades automatically reduced by 50%. 

Syntax (or any other) errors: If you program does not run (i.e., if it returns any kind of error 
message) your grade will be automatically reduced by 50%. 



